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INTRODUCTION 

The  following  five  algorithms  formulate  the  perturbing  potentials  associated 
with  the  four  Earth  tides^ * 3,4 , 5  terms  of  linear  combinations  of  the  eigen¬ 

functions  Unm  and  Vnm  (solid  harmonics)  used  to  express  the  main  gravitational 
field  in  OMNIS. ^  They  make  it  possible  to  calculate  the  gravitational  effect  on 
the  satellite  orbit  of  the  mass  redistribution  caused  by  the  tides  by  simply 
adding  the  generally  time-dependent  expansion  coefficients  of  the  expressions  for 
the  tidal  perturbing  potentials  as  increments  to  the  expansion  coefficients  Cnm 
and  Snm  of  the  main  gravitational  potential  field.  Also,  by  processing  the  tidal 
increments  to  the  main  field  coefficients  along  with  the  latter,  the  computation 
of  the  second  derivatives  with  respect  to  the  field-point  coordinates  necessary  to 
include  the  contribution  of  the  Earth  tides  to  the  variational  equations  occurring 
in  connection  with  the  orbit  determination  becomes  an  entirely  automatic  part  of 
the  contribution  by  the  main  gravitational  field  parameters  to  the  variational 
equations . 

The  first  and  second  algorithms  specify  the  lunar  atmospheric  tide  and  the 
air  tides  produced  by  the  combined  gravitational  and  thermal  effects  of  the  Sun, 
respectively.  These  algorithms  had  originally  been  formulated^ * ^ * ®  directly  as 
perturbing  accelerations.  They  were,  for  the  purposes  of  OMNIS,  rewritten  in 
terms  of  the  required  potential  expansion  coefficient  increments.  Only  minor 
changes  were  needed  to  adapt  the  CELEST  Extended  Ocean  Tide^  to  OMNIS  as  its  third 
Earth  tide  routine.  The  desired  expansion  coefficient  increments  were  already 
part  of  the  CELEST  version.  The  fourth  OMNIS  tide  algorithm  accounts  for  the 
solid  Earth  tide  as  specified  by  the  MERIT  campaign  document.-^  This  particular 
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tide  formulation  was  readily  expressed  in  the  form  necessary  for  OMNIS.  And, 
finally,  the  fifth  routine  is  a  reformulation  of  the  familiar,  elementary,  solid 
Earth  tide  that  was  part  of  earlier  orbit  integration  routines  contained  in 
computer  programs  such  as  CELEST,^  TERRA, and  ASTRO. ^  It  is  included  here 
because  some  users  of  OMNIS  may  occasionally  want  to  activate  it. 

DATA  ASSOCIATED  WITH  TIME  AND  THE  FUNDAMENTAL  ASTRONOMICAL  ARGUMENTS 

INPUT  DATA 

Y  year  (four  digits) 

D  day  number  (number  of  day  in  year  Y) 
t*  Universal  Time  in  seconds  (s) 

ALGORITHM 

Using  double  precision  throughout,  evaluate  the  time  instant  T  of  the  time 
line  under  consideration,  expressed  in  Julian  Centuries  of  36525  days  of  86400  s 
of  Dynamical  Time  elapsed  since  Epoch  J2000.0  (=  JD  2451545.0  =  2000  January  1.5) 

M  =  1  (101) 

Nj.  =  367  *  Y  -  7  *  (Y  +  (M  +  9)/12)/4 
-  3  *  ((Y  +  (M  -  9 )/7 )/100  +  1 )/4 

+  275  *  M/9  +  1721029  +  D  (102) 

MJDO  =  Nz  -  0.5  -  2400000 

=  Modified  Julian  Date  at  beginning  of  day  D,  UT 

T  _  MJDO  -  51545 
0  36525 


(103) 

(104) 
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T  =  Tq  + 


86400  *  36525 


(105) 


Note  that  the  equation  for  N^,  is  a  FORTRAN  expression.  When  evaluating  it  outside 
of  a  FORTRAN  computer  program,  the  rules  of  FORTRAN  integer  arithmetic  must  be 
observed.  This  expression  computes  the  Julian  day  number  corresponding  to  a  time 
instant  43200  s  into  (noon  UT  of)  the  day  D,  during  month  M  and  year  Y.  It  is 
valid  for  all  A.D.  years.  Y,  M,  and  D  indicate  the  Gregorian  calendar  date. 
Evaluate  also 


l  =  134.96298139  +  477198.867398056  T 

+  0.008697222  T2  +  0.000017778  T3 

l'  =  357.527723333  +  35999.050340000  T 

-  0.000160278  T2  -  0.000003333  T3 

F  =  93.271910278  +  483202.017538056  T 

-  0.003682500  T2  +  0.000003056  T3 

D  =  297.850363056  +  445267.111480000  T 

-  0.00191417  T2  +  0.00000528  T3 

fl  =  125.0445222  -  1934.13626083  T 

+  0.00207083  T2  +  0.00000222  T3 


(106) 


(107) 


(108) 


(109) 


(110) 


All  five  quantities  result  in  decimal  fractions  of  degrees.  Further  evaluate  the 
lunar  time  angle,  in  degrees,  reckoned  from  lower  transit 


x  =  0g  +  180°  -  s 


the  Moon's  mean  longitude,  in  degrees 


s  =  f  +  a 


The  Sun's  mean  longitude,  in  degrees 


(111) 


(112) 
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the  longitude,  in  degrees,  of  the  Moon's  mean  perigee 


p  =  s  -  l 


(114) 


the  negative  longitude,  in  degrees,  of  the  Moon's  mean  node 


n  =  -  n 


(115) 


and  the  longitude,  in  degrees,  of  the  Sun's  mean  perigee 


PX  =  s  -  D  -  l’ 


(116) 


where 


9  is  the  Greenwich  mean  sidereal  time,  in  degrees. 

o 


9  =  100.4606184  +  36000.7700537  TQ 

8  +  0.000387933  Tq 


1.002737909 

+  -  •  ]S  ♦ 

3600 


(117) 


Reduce  the  quantities  calculated  by  Equations  106  through  117  Modulo  360. 


SPECIFICATION  OF  EIGENFUNCTIONS 


The  following  equations  are  not  part  of  the  algorithm.  They  are  statements 


intended  to  clarify  what  precisely  is  meant  whenever  the  term  "eigenfunction" 


appears  in  the  present  coding  request, 


L  =  u  - r  P  (sin  9)  cos  mX 

n  n+1  n 

r 


(201) 


V  -  - -  P  (sin  v)  sin  mA 

n  n+ 1  n 

r 


(202) 


VS* 


smews* 
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where 

u  =  Earth's  gravitational  parameter  (=  gravitational  constant  G 
times  Earth’s  mass) 

R  =  semimajor  axis  of  reference  ellipsoid  (Earth's  "radius") 
r  =  geocentric  distance  of  field  point 
0  =  field  point  geocentric  latitude 

X  =  field  point  longitude  in  Earth-fixed  frame 

p™(x)  =  Legendre  Polynomials  (normalization  as  stated  in  "JAHNKE-EMDE" ) 


The  relationship  of  these  eigenfunctions  with  the  gravitational  expansion 
coefficients  Cn>fn  and  Snjm  and  with  the  Earth's  main  gravitational  potential,  U, 
is 


U  =  Z  (c 

n,m 


n,m 


ura  +  s  vm) 

n  n  ,m  n 


(203) 


The  following  computer  program  segments  correspond  to  the  individual  tides 
and  produce  increments  ACnjin  and  ASn<m  to  the  CntIT,  and  Snjin  ,  These  increments 
reflect  the  effect  of  the  particular  Earth  tide  on  the  satellite  motion  and  on  the 
variational  equations. 

t_ue  that  the  above  eigenfunctions,  potential  expansion  coefficients,  and 
increments  to  the  latter  expansion  coefficients  are  those  quantities  that  are 
normally  referred  to  as  "unnormalized"  eigenfunctions  and  expansion  coefficients 
to  be  distinguished  from  the  corresponding  "normalized"  quantities.*  Note  that 
the  MERIT  Campaign  Document  uses  normalized  quantities  for  its  version  of  the 
solid  Earth  tide. 


*If  an  elaboration  of  the  terms 
NSWC  TR-2223  A  Geoid  Height  and 
Parameter  Solutions.  W.  Groeger, 


"normalized"  and  "unnormalized"  is  needed, 
Gravity  Anomaly  Computer  Program  For  Geodet 
among  other  references. 


see 

ic 
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As  far  as  units  and  dimensions  are  concerned,  it  is  assumed  that  for  the 
purposes  of  formulating  the  Earth’s  main  gravitational  field,  4 ,  R,  and  r  will  be 
expressed  in  terms  of  km  and  s  (4  =  398600,5  km^  s-^ ) .  Accordingly,  all  C  and  S 
increments  occurring  below  will  be  equipped  with  numerical  factors  that  ensure  the 
values  calculated  from  these  expressions  will  be  compatible  with  the  expansion 
coefficients  of  the  main  gravitational  potential  field. 


CONTRIBUTIONS  TO  THE  GRAVITATIONAL  POTENTIAL  EXPANSION 
COEFFICIENTS  BY  THE  LUNAR  ATMOSPHERIC  TIDE 


INPUT  DATA 


A2  =  a  constant  associated  with  the  physics  of  the  lunar  air 

tide.  Suggested  value  A2  =  0.564  kg  m 

R  =  any  reasonable  value  for  the  mean  "radius”  of  the  Earth 

(semimajor  axis  of  the  earth's  reference  ellipsoid), 
in  km.  Default  value  R  =  6378.140  km 

G  =  gravitational  constant  in  m,  kg,  and  s. 

Suggested  value  G  =  6.6732E-11  m^  kg~^  s-^ 


4  =  Earth's  gravitational  parameter. 

Suggested  value  4  =  398600.5  km^  s-^ 

Y  =  year 


D  =  day 

t*  =  Universal  Time  in  s 


s(T)  and  h(T)  =  fundamental  astronomical  parameters 


ALGORITHM 


a 


=  A2  G  R 


(301) 


l  =  a/48 


(302) 


t** 


360 

86400 


t* 


(303) 
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v  =  s(T)  -  h(T)  (304) 


r  =  2(t**  -  v)  -  15°  (305) 


The  increments  to  the  gravitational  potential  expansion  coefficients,  due  to 


the  lunar  air  tide,  are 


AC2.2  = 

+  ** 

V 

10"3 

cos 

r 

(306) 

AS2,2  = 

_  aR 
M 

10'3 

sin 

r 

(307) 

AC4  f 2  = 

_  bR 
y 

10~3 

cos 

r 

(308) 

AS4 1 2  = 

+  M 

y 

10-3 

sin 

r 

(309) 

CONTRIBUTIONS  TO  THE  GRAVITATIONAL  POTENTIAL  EXPANSION 
COEFFICIENTS  BY  THE  SOLAR  ATMOSPHERIC  TIDE 

INPUT  DATA 

=  a  constant  associated  with  the  physics  of  the  solar  air  tide. 

Suggested  value  Bj  =  6  kg  m”2 

B2  =  a  constant  associated  with  the  physics  of  the  solar  air  tide. 

Suggested  value  B2  =  11.9  kg 

R  =  semimajor  axis  of  Earth's  reference  ellipsoid,  in  kilometers. 

Default  value  R  =  6378.140  km 

G  =  gravitational  constant  in  m,  kg,  and  s. 

Suggested  value  G  =  6.6732E-11  m^  kg“l  s-^ 

U  =  Earth's  gravitational  parameter.  Suggested  value  y  =  398600.5  km^  s-^ 


t*  =  Universal  Time  in  s 
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ALGORITHM 


al  =  B1 


8ttGR 

105 


a2  *  b2 


5tt2GR 

64 


!  a3  =  b2 

i 


5tt2GR 

3072 


=  a2/48 


The  increments  to  the  gravitational  potential  expansion  coefficients, 
to  the  solar  air  tide,  are 


alR  -3 

in  J 


360 


AC3,1  -  10  cos  (g^  t*  -  78°) 

a  R 

ss3.i  *  +  -r  10'3  sln  ‘sUbo  -  78°> 


a2R  -3 

•-  1  A  J 


ac2,2  =  +  10  cos  2  ( 


32R  -3 

in  J 


360 

86400 


360 


t*  -  146°) 


aS2,2  -  10  sin  2  <8^50  C*  "  U6°> 

a  R 

“C/,,2  •  -  -jr  10-3  cos  2  (gffgg  t*  -  U6') 


a3R  _3 

1  n  J 


360 


^4,2  =  +  ~ -  10  sin  2  t*  -  146°) 


CONTRIBUTIONS  TO  THE  GRAVITATIONAL  POTENTIAL  EXPANSION 
COEFFICIENTS  BY  THE  OCEAN  TIDE 


INPUT  DATA 


Concerned  are  the  ocean  tide  modes  M2,  S2,  N2,  K2,  Kj,  O3,  Pj ,  Q3 ,  Mf 
and  Ssa.  For  each  of  these  tides,  there  exists  a  computer  tape  that  lists 
values  of  the  tidal  amplitude  r  3  j  and  the  tidal  phase  63 j . 


(401) 

(402) 

(403) 

due 

(404) 

(405) 

(406) 

(407) 

(408) 

(409) 


« 

the  NP 


NP  Number  of  grid  points.  This  may  be  expected  to  be  a  five-digit  integer. 

Cij  Tidal  amplitude  on  area  element  ij,  If  available  in  meters,  use 
directly.  Otherwise,  convert  to  meters. 

Tidal  phase  angle  on  area  element  ij,  in  degrees. 

Custodian  of  the  tapes  just  mentioned  is  Ling  Szeto  of  the  Physical  Sciences 
Software  Branch.  Also  required  are 

R  "Radius  of  the  Earth,"  i.e.,  the  semimajor  axis  of  a  suitable 
reference  ellipsoid,  in  km.  Trial  value  R  =  6378. 145  km. 

Square  of  eccentricity  of  reference  ellipsoid. 

Trial  value  =  0.00669342,  In  case  it  is  desired  to  start  from  the 
ellipsoid  flattening,  f,  find  from  =  (2  -  f)f. 

P  Gravitational  constant  of  the  Earth,  in  km  and  s. 

Trial  value  p  =  398600.5  km^  s-^ 

G  Newton's  constant,  in  km,  kg,  and  s. 

Trial  value  G  =  6.6732E-20  km^  kg-^  s-^. 

P  Density  of  water,  in  kg  and  km. 

Trial  value  0  =  l.E  +  12  kg  km~^. 

P*  Density  of  the  ocean  bottom,  in  kg  and  km. 

Trial  value  p*  =  3.E  +  12  kg  km-^. 

o  Rate  of  mean  longitude  of  Moon. 

o  =  —  1 . 40519E-04  deg  s_1 . 

TT 

NMAX  Limit  on  the  number  of  terms  in  the  expansion  of  the  tide  potential. 

The  following  parameters  Y,D,  and  t*  specify  the  time  instant  at  which  the  per¬ 
turbing  acceleration  due  to  the  tide  is  to  be  found.  Normally,  this  is  the  time 
instant  associated  with  the  current  time  line  of  the  orbit  integration. 

Y  Number  of  the  calendar  year 

D  Number  of  the  day  within  the  year 

t*  Mean  solar  time  at  Greenwich  (GMT,  UTC),  in  s 
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ALGORITHM 


Computer  Routine  for  the  Expansion  Coefficients  of  the  Tide  Potential 

Ct  fi 

The  time-independent  constituents  of  the  expansion  coefficients  (  Fnm,  pFnm, 

a  8 

Hnm,  Hnn))  of  the  tide  potential  will  now  be  calculated  from  the  amplitude  and 
phase  angle  computer  tapes,  separately  for  each  of  the  11  ocean  tide  modes. 

Once  established,  the  coefficients  Fnm,  PFnm,  Hnm,  Hnm  will  remain 
unchanged  as  long  as  the  amplitude  and  phase  tape  for  the  particular  tide  mode 
remains  valid.  Thus,  the  computer  routine  under  discussion  will  be  exercised 
quite  infrequently,  namely  just  once  each  time  a  new  edition  of  the  just-mentioned 
amplitude  and  phase  tape  becomes  available.  That  is  expected  to  occur  at  inter¬ 
vals  of  several  years,  during  which  the  present  routine  will  remain  dormant. 

For  each  ocean  tide  mode  execute  the  following  algorithm.  For  the  geometry 
associated  with  the  index  ij,  see  Figures  1,  2,  and  3. 


C4j  =  10“3(P  -  0.0667P*)GASij£ijCOs6ij  (501) 

Sij  =  10~3(P  -  0.0667p*)GASijCijSin6ij  (502) 

2 

ASij  =  (l8o)  R  Sin  (l80  (503) 


=  area  of  the  nonpolar  surface  area  element  for 
j  =  2, 3,4,  •  •  ’ ,  jj^x  <  180 


(ASijJpQLAR 


2 


3 


R 


2 


(504) 


=  area  of  the  polar  surface  area  element  neighboring  the 
North  Pole.  (Note  that  the  southern  polar  surface  area 
element  is  not  needed  because  the  computer  tape  mentioned 
on  page  9  contains  tide  data  for  latitudes  above  60  deg 
South  only. ) 
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A  =  li-1) 


A  =  A.  =  i 


FIGURE  2.  POSITION  OF  POINT  MASS  AT  THE  GEOMETRICAL  CENTER 
OF  THE  SURFACE  AREA  ELEMENT  ASSOCIATED  WITH  ij 
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v,  -  INDICATES  THE  POINT  WHERE  GREENWICH  MERIDIAN  AND 
EQUATOR  INTERSECT 
v2  *  LOCATED  IN  THE  EQUATORIAL  PLANE 


FIGURE  3.  POSITION  OF  THE  POINT  MASS  M ij  IN  THE 
EARTH-FIXED  CARTESIAN  COORDINATE  FRAME 


\  ~  155  (J  "  b  ln  rad 


Bo  (1  ’  1>  ln  rad 


(i  -  y  si°2  9ij) 


To  each  index  tj  assign  now  a  counting  number  v,  1  <  v  <  NP. 


'  <2  -  Sm> 

(n  - 

m)  ! 

1 

NP 

V  o2n+1 

(n  +  m)  ! 

D211 

wR 

l  v 

v=  1 

=  (2  -  6°) 
m 

(n  - 

m)  ! 

1 

fp  2n+l 
v=  1 

(n  + 

m)  ! 

„2n 

ER 

(n  -  m)  ! 

1 

NP 

£ 

V=1 

2n+l 

(n  +  m)  ! 

u  R2n 

pv 

ot  n 
v  nm 

(n  -  m)  ! 

1 

NP 

£ 

V=1 

2n+l 

2  hV 

v  nm 

(n  +  m)  ! 

2n 

uR 

Dv 

(505) 


(506) 


(507) 


(508) 


(509) 


(510) 


(511) 


\V 


i-V*- 

1-  9 


rJf  V 


AiI  »%  % 
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rk  1  f  l  -  k 

&i  0  if  l  i  k 

v  v 

Calculate  the  fnm  and  hnm  as  follows,  noting  that 


Pv  = 


(512) 


(513) 


and 


hno  =  0  for  all  values  of  n 


(514) 


and  that 


hnm  =  fnm  =  0  for  anX  n  <  m‘ 


(515) 


Obtain,  separately  for  each  v,  the  required  fpm  and  h^Jm  from  the  following 
recurrence  relations: 

To  advance  in  n,  evaluate 
P 


n+l,m  n  -  m 
P 

v  _ v_ 

^n+l,m 


T-r  [  (2n  +  1)  sin  9  f  -  (n  +  m)p  f  ] 
+  1  v  n,m  v  n,m 


n  -  m 


— — 7  [(2n  +  1)  sin  0  h  -  (n  +  m)p  h  ] 

+  1  v  n,m  v  n-1 ,m 


(516) 

(517) 


To  advance  in  m,  use 


fv  ,  ,  =  (2n  +  l)p  [cos  9  cos  A  fV  -  cos  0  sin  A  hV  ] 

n+l,n+l  v  v  v  n,n  v  v  n,n 

hV  ,  ,  =  (2n  +  l)p  [cos  8  cos  A  hV  +  cos  6  sin  A  fV  ] 

n+l,n+l  *v  v  v  n,n  v  v  n.n 


(518) 

(519) 


Start  the  recurrences  from 


fV  =  — 

°,°  pv 


sin  9 


£1.0  ■  R 


h0,0  '  hl,0  ■  0 


(520) 

(521) 

(522) 


and  terminate  the  procedure  at  n  =  NMAX. 
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Store  the  resulting  “Fnn,,  ^Fnin,  “H,,,,,,  and  ®Hnm.  They  will  be  programmed  as 
constant  parameters  into  the  computer  routine  for  the  perturbing  potential  and  are 
expected  to  serve  for  all  subsequent  computer  runs  of  that  routine,  until  updated. 
Note  that  the  aFnm,  ^Fnm,  aHnm,  and  ®Hnm  are  dimensionless  quantities. 


Calculation  of  the  Astronomical  Arguments  for  the  Various  Tide  Components 


Evaluate  now,  using  double  precision  throughout,  the  Moon's  longitude,  s,  the 
Sun's  mean  longitude,  h,  and  the  longitude  of  the  Moon's  mean  perigee,  p,  at  the 
beginning  of  the  day  UT  in  question 


s0  =  s(Tq) 
h0  =  h(T0) 

PO  =  P<T0) 

For  each  tide  mode,  find  now  the  Astronomical  Argument  X  from 


(523) 


(524) 


(525) 


Tide  Mode 


X  in  deg 
2 (hQ  ~  s0) 


2h0  -  3S0  +  PO 


h0  +  90 

h0  -  2S0  -  90 

-ho  -  90 

h0  -  3S0  +  PO  -  90 


so  -  PO 


(526) 


(527) 


(528) 


(529) 


(530) 


(531  ) 


(532) 


(533) 


(534) 


(535) 


(536) 


%y, 

w. 

iVV 


•  *■< 

«»*; 


m 

1 


!» -.f 

*4 

ipS 

h'.‘ 


I 


iW*. 


Si* 

Ss 

& 
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Note  that  the  day  number  D,  and  thus  X,  must  be  updated  whenever  the  time 
argument  t*  runs  into  the  next  day. 

Computer  Routine  for  the  Increments  to  the  Gravitational  Potential  Expansion 
Coefficients 


For  each  tide  mode,  separately,  use  the  above  set  of  time-independent  tide  M 
potential  coefficients  and  the  astronomical  tide  argument  and  calculate  the  time-  in 
independent  expansion  coefficients  for  the  tide  potential  (note  n,m  <  NMAX).  The  |1 
latter  are  the  increments  to  the  gravitational  potential  expansion  coefficients  I 
due  to  the  particular  ocean  tide  mode  9 


^Cnni 

=  aFnmcos(at*  +  X) 

+  8Fnmsin<ot* 

+ 

X) 

(537) 

ASnm 

=  aHnmcos(at*  +  x) 

+  3Hnmsin(ot* 

+ 

X) 

(538) 

The  argument  (at*  +  X)  should  be  evaluated  in  double  precision.  The  double  pre¬ 
cision  sin  and  cos  functions  should  be  applied.  Then  calculate  -Fnm  and  AHnm  and 
subsequently  revert  to  single  precision. 

Note  also  that  AFnm  and  AHnra  are  linear  functions  of  the  two  trigonometric 
terms.  To  evaluate  these  trigonometric  terms,  let  t^  and  t*+j  be  the  values  of 
Universal  time  for  which  subsequent  integration  steps  are  to  be  performed.  To 
save  computer  time,  update  the  trigonometric  time  factors  as  follows 

cos(ot*+i  +  x)  =  cos[(ot*  +  x)  +  oAt] 

=  cos(ot*  +  x)cos  aAt  -  sin(at*  +  x)sin  a At  (539) 

sin(at*+i  +  x)  =  sinfat*  +  x)  +  oAtJ 

^  & 

=  sin(  ;ti  +  x)cos  oAt  +  cos(atj  +  x)sin  oAt  (540) 


9 


4 


% 
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CONTRIBUTION  TO  THE  GRAVITATIONAL  POTENTIAL  EXPANSION 
COEFFICIENTS  BY  THE  SOLID  EARTH  TIDE 


INPUT  DATA 


Y  year 

D  number  of  day  in  year 

t*  Universal  Time  in  s 

k2  nominal  second  degree  Love  number.  Default  value  k£  =  0.3. 

R  semimajor  axis  of  reference  ellipsoid  in  km. 

Default  value  R  =  6378.140  km 

P  Earth’s  gravitational  parameter,  p  =  398600.5  km-*  s-^ 

p2  Moon's  gravitational  parameter.  P2  =  4916.816  km^  s-^ 

P^  Sun’s  gravitational  parameter.  P3  =  132.712E+09  km^  s-^ 

x*  \  the  inertial,  Cartesian  coordinates  of  the  Moon,  in  km, 

I  for  the  time  instant  (Y,  D,  t*),  from  the  JPL  Export 
y£  \  Ephemeris.  The  reference  frame  is  referred  to  equator 
I  and  equinox  corresponding  either  to  Epoch  J2000.0  or  to 
z*  I  Epoch  1950.0. 

x*  I  the  inertial,  Cartesian  coordinates  of  the  Sun,  in  km,  for 

*  I  the  time  instant  (Y,  D,  t*),  from  the  JPL  Export  Ephemeris. 

^3  /  The  reference  frame  is  referred  to  equator  and  equinox 

*  1  corresponding  either  to  Epoch  J2000.0  or  to  Epoch  1950.0. 


ALGORITHM 


The  Earth-Fixed  Moon  and  Sun  Coordinates 

Calculate  now  the  Earth-fixed  coordinates  of  the  Moon  (j  =  2)  and  of  the 
Sun  (j  =  3),  Xj,  yj  and  z j ,  at  the  time  instant  (Y,  D,  t*)  from  the  corresponding 
inertial  coordinates 


(ABCD) 


Y.D.t* 


(601) 
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where  (ABCD)ytn<t*  is  the  transformation  matrix  that  manages  the  transition  from 
the  (inertial)  ephemeris  frame  of  reference  to  the  Earth-fixed  reference  frame 
associated  vitji  the  time  instant  (Y,  D,  t*).  This  transformation  may  be  assumed 
to  be  a  readily  available  subroutine  of  OMNIS. 

Further  calculate,  for  j  =  2  and  j  =  3, 


ri  -  M 

12.  2  2 
'  X  ,  +  y  ,  +  z  . 

J  J  J 

=  +  1 

/(X*)2  ♦  <y*>2  +  O*)2 

(602) 

Pj  =  +  1 

/  2  .  2 
/x.  +  y. 

'  J  J 

(603) 

COS  <j>j  = 

pj/rj 

(604) 

sin  <j>j  = 

zj/rj 

(605) 

sin  2 4> j 

=  2  sin  4» j  cos  4>j 

(606) 

cos  Xj  = 

xJ/pj 

(607) 

sin  Xj  = 

y  j/p  j 

(608) 

cos  2Xj 

=  cos2Xj  -  sin2X^ 

(609) 

sin  2X j 

=  2  sin  X  .  cos  X  . 

J  J 

(610) 

where 


rj  =  distance  from  the  geocenter  to  the  Moon/Sun  in  km 
=  geocentric  latitude  of  the  Moon/Sun 
> j  =  longitude  (East)  of  the  Moon/Sun 


The  Long-Period  Tide  Correction 

Compute  the  increments  for  the  normalized  geopotential  expansion  coefficient 


C20  as  follows  . 
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,3  3 


AC 


20 


^k2  v  £  ,_\3  P20(sin  'V 


j-2  (Tj ) 


20' 


P2Q(x)  -  |  (3  x2  -  1) 


(611) 

(612) 


The  term  <(AC2o)>is  *  constant  reflecting  the  "Permanent  Tide" 

<AC20>  =  -  1.39119  •  lO-8  k2  (613) 

It  must  be  deleted  when  using  C2o  values  derived  by  a  method  that  already  accounts 
for  the  effect  of  the  Permanent  Tide. 


The  Diurnal  Tide  Corrections 


Calculate  now  the  increments  to  the  normalized  geopotential  expansion  coeffi¬ 
cients  C21  and  S21  (n  =  2,  m  =  1,  n  +  m  =  odd) 

,3  3 


AC. 


-IJI  k  ± 

!  1  3  T5  k2  u  A 


P  j 


3=2  (r.) 


—  P21(sin  4>j)cos  A  +  AAC21 


P21  (sin  4>j)  =  3  1  -  sin2  <p ^  sin  <p j  =  sin  2  4>j 

AAC2i  =  Y'  (A  6k  H  )sin  0 
mss  s 


16.4 

.  10-12 

sin 

(  T  -  S) 

49.6 

•  10~12 

sin 

(x  +  s- 

9.4  • 

10~12 

sin 

(  T  +  S  -  N 

507.4 

•  10-12  sin  (t  +  s) 

73.5 

•  10"12 

sin 

(  I  +  S  + 

15.2 

.  10-12 

sin 

(  T  +  S  + 

16.4 

•  io-12 

sin 

(0  -  2F 

g 

49.6 

•  10~12 

sin 

(0g  -  2F 

9.4  • 

10-12 

sin 

(0g  +  2) 

507.4 

•  10  I2  sin  6g 

73.5 

•  10-12 

sin 

(9g  -  U) 

15.2 

•  10~12 

sin 

<°8  *  r> 
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(614) 

(615) 


(6 16A) 


(616B) 
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AS21 


7^  P21<sl”  *j)sin 

(r.) 


X  . 
J 


+  AAS 


21 


AAS2l  =  16.4  •  lO"12  cos  (0g  -  2F  -  2^1 ) 

+  49.6  •  lO"12  cos  (0g  -  2F  -  20  +  2D) 

+  9.4  •  lO-12  cos  (0g  +  Q) 

-  507.4  •  10”12  cos  0g 

-  73.5  •  10-12  cos  (0g  -  a) 

+  15.2  •  10"12  cos  (0g  +  l') 


(617) 


(618) 


The  Semidiurnal  Tide  Corrections 


Calculate  finally  the  increments  to  the  normalized  geopotential  expansion 
coefficients  C22  and  ^22  (n  =  2,  m  =  2,  n  +  m  =  even) 


ac22  =  -iWr  k2  t  4  77b  P22<sin  Vcos  “j  + 

J  '•rj ' 

P22  (sin  <fij)  =  3  cos2  <(> j 

AAC22  =  39.5  •  10"12  cos  2(0g  -  F  -  a) 

+  18.4  •  10~12  cos  2(0g  -  F  -  a  +  D) 

“22  =  T?  Vt  k2  T  77^  P22(Sln  V”1"  2>j  + 

AAS22  =  -  39.5  *  10-12  sin  2(9g  -  F  -  u) 

-  18.4  •  10-12  sin  2  ( 6g  -  F  -  -l  +  D) 


(619) 

(620) 

(621) 

(622) 

(623) 


The  increments  to  the  gravitational  potential  expansion  coefficients,  due  to  the 
solid  Earth  tide,  are 


ACn,m  =  AC 


n,m  f 


m)  !  (2n  +  1)  (2  -  6  ) 
_ m 

(n  +  m)  ! 


(624) 
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CONTRIBUTION  TO  THE  GRAVITATIONAL  POTENTIAL  EXPANSION  COEFFICIENTS  BY  THE 
SOLID  EARTH  TIDE  (TERRA/CELEST  FORMULATION— BACKUP  FORMULATION  FOR  OMNIS) 


INPUT  DATA 


Y 

D 

t* 

k2 

R 


year 

number  of  day  in  year 
Universal  Time  in  s 

nominal  second  degree  Love  number.  Default  value  k2  =  0.3. 


semimajor  axis  of  reference  ellipsoid  in  km. 
Default  value  R  =  6378.140  km. 


398600.5  km3  s'2 , 


Earth's  gravitational  parameter.  V 
Moon’s  gravitational  parameter.  t>2  =  4916.816  km3  s-2 . 
Sun's  gravitational  parameter.  U3  =  132.712E+09  km3  s~2 , 


the  inertial,  Cartesian  coordinates  of  the  Moon,  in  km, 
for  the  time  instant  (Y,  D,  t*),  from  the  JPL  Export 
Ephemeris.  The  reference  frame  is  referred  to  equator 
and  equinox  corresponding  either  to  Epoch  J2000.0  or  to 
Epoch  1950.0, 


the  inertial,  Cartesian  coordinates  of  the  Sun,  in  km,  for 
the  time  instant  (Y,  D,  t*),  from  the  JPL  Export  Ephemeris. 
The  reference  frame  is  referred  to  equator  and  equinox 

corresponding  either  to  Epoch  J2000.0  or  to  Epoch  1950.0. 
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ALGORITHM 

The  Earth-Fixed  Moon  and  Sun  Coordinates 

Calculate  now  the  Earth-fixed  coordinates  of  the  Moon  (j  =  2)  and  of  the 
Sun  (j  =  3),  xj,  yj  and  zj,  at  the  time  instant  (Y,  D,  t*)  from  the  corresponding 
inertial  coordinates 


where  (ABCD)y^D,t*  *s  the  transformation  matrix  that  manages  the  transition  from 
the  (inertial)  ephemeris  frame  of  reference  to  the  Earth-fixed  reference  frame 
associated  with  the  time  instant  (Y,  D,  t*).  This  transformation  may  be  assumed 
to  be  a  readily  available  subroutine  of  OMNIS. 

Further  calculate,  for  j  -  2  and  j  =  3 


V2  ,  2  ,  2 

x  .  +  y  .  +  z  . 

J  J  J 


=  +  y (x*)2  +  (Y*)2  +  (z*)2 

(702) 

/  2  ,  2 
j  "  J  3 

(703) 

cos  4> j  =  Pj/rj 

(704) 

sin  4j  =  z j/r j 

(705) 

sin  24j  =  2  sin  <J>j  cos  <t>j 

(706) 

COS  ‘j  =  Xj/Pj 

(707) 

sin  =  yj/Pj 

(708) 

cos  2‘j  =  cos^A  -  sin^A 

(709) 

sin  2j  =  2  sin  A  cos  A 

(710) 
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where 

r j  *  distance  from  the  geocenter  to  the  Moon/Sun  in  km 
♦j  a  geocentric  latitude  of  the  Moon/Sun 
Aj  ■  longitude  (East)  of  the  Moon/Sun 

The  Increments  to  the  Potential  Expansion  Coefficients 


The  increments  to  the  gravitational  potential  expansion  coefficients,  due  to 
the  solid  Earth  tide  (TERRA/CELEST  Version),  are 


^2,0  = 


^2,1 


AS2.1  = 


AC2.2  * 


AS2,2 


3  y.  r3 

k2  E  -?■—  P2(sln  ♦j) 

J-2 

k2  E  "S  t  )c»»  X 

J-2  r. 

3  .  v .  r3  . 

k2  I  1  P2(sin  4>  )sin  A 

i=2  r.  J  J 

J  i 


3  I  P  •  r3  o 

k2  In  f  n  p2(sln  Vcos  2Xj 

k2  I [72  f  -S  P2(sl"  *j)sln  2Aj 

J'2  rj 


(711) 


(712) 


(713) 


(714) 


(715) 


where 


?2  (sin  )  -  yO  sin^  $  _  d 


?2  (sin  $j)  H  3  -y  1  -  sin^  sin  4>j  =  y  sin  24>j 
2 

?2  (sin  4>j)  =  3  cos^  ^ 


(716) 


(717) 


(718) 


Note  that  the  present  algorithm  is  a  special  case  of  the  algorithm  for  the  con¬ 
tribution  to  the  gravitational  expansion  coefficients  by  the  solid  Earth  tides 
on  page  16  (Wahr  theory  terms  deleted). 


V  V.  V 
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